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Abstract 

We investigate the prompt photon photoproduction at HERA within the framework of h?- 
factorization QCD approach. Our consideration is based on the off-shell matrix elements for 
the underlying partonic subprocesses. The unintegrated parton densities in a proton and in a 
photon are determined using the Kimber-Martin-Ryskin (KMR) prescription. Additionally, 
we use the CCFM-evolved unintegrated gluon as well as valence and sea quark distributions 
in a proton. A conservative error analisys is performed. Both inclusive and associated with 
the hadronic jet production rates are investigated. The theoretical results are compared with 
the recent experimental data taken by the HI and ZEUS collaborations. We study also the 
specific kinematical properties of the photon-jet system which are strongly sensitive to the 
transverse momentum of incoming partons. Using the KMR scheme, the contribution from 
the quarks emerging from the earlier steps of the parton evolution is estimated and found 
to be of 15 - 20% approximately. 

PACS number(s): 12.38.-t, 13.85.-t 



1 Introduction 

The prompt photon production in ep collisions at HERA is subject of intense studies [1- 
6]. The theoretical and experimental investigations of such processes have provided a direct 
probe of the hard subprocess dynamics, since produced photons are largely insensitive to the 
effects of final-state hadronization. Usually photons are called "prompt" if they are coupled 
to the interacting quarks. From the theoretical point, these photons in ep collisions can be 
produced via direct — » -yq and resolved production mechanisms. In resolved events, the 
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photon emitted by the electron fluctuate into a hadronic state and a gluon and/or a quark of 
this hadronic fluctuation takes part in the hard interactions. Prompt photon measurements 
can be used also to constrain the parton densities in the proton and in the photon. 

Recently the HI and ZEUS collaborations have reported data [2-6] on inclusive and 
associated (with the hadronic jet) prompt photon production at HERA. However, next- 
to-leading order (NLO) collinear pQCD calculations [7, 8] are 30 — 40% below these data, 
especially in rear pseudo-rapidity (electron direction) region. It was demonstrated [2-5] that 
the observed disagreement is difficult to explain with conventional theoretical uncertainties 
connected with scale dependence and parametrizations of the parton densities. The origin of 
the disagreement has been ascribed to the effect of initial-state soft-gluon radiation. It was 
shown [3] that observed discrepancy can be reduced by introducing some additional intrinsic 
transverse momentum kr of the incoming partons. The ZEUS fit to the data gave a kx value 
of about 1.7 GeV [3]. A similar situation is observed also at Tevatron energies: in order to 
describe the measured transverse momentum distributions of the photon the Gaussian-like 
kr spectrum with an average value of kx ~ 3 GeV was introduced [9, 10]. Of course, such 
large partonic kx must have a significant perturbative QCD component. 

The transverse momentum of incoming partons naturally occurs in the framework of 
A; T -factorization approach of QCD [11]. In this approach, the transverse momentum k T is 
generated perturbatively in the course of non-collinear parton evolution via the corresponding 
(usually Balitsky-Fadin-Kuraev-Lipatov (BFKL) [12] or Ciafaloni-Catani-Fiorani-Marchesini 
(CCFM) [13]) evolution equations. A detailed description of the /cy-factorization can be 
found, for example, in reviews [14-16]. As it was demonstrated in the ZEUS paper [4] and 
in the recent experimental study [6] performed by the HI collaboration, the /^-factorization 
predictions [17] for prompt photon photoproduction at HERA are in better agreement with 
the data than the published results of the collinear NLO pQCD calculations [7, 8]. 

An important component of the first calculations [17] in the framework of /^-factorization 
approach was the unintegrated quark distributions f q (x, k. T , /z 2 ) in a proton. These quantities 
are poorly known since there are theoretical difficulties in obtaining the quark distributions 
directly from CCFM equation (see also [14-16] and references therein for more informa- 
tion). At present, the unintegrated quark densities are most often used in the framework of 
KMR [18] approximation only. As a result, the dependence of the /c^-factorization predic- 
tions [17] on the non-collinear evolution scheme has not been investigated. This dependence 
in general can be significant and it is a special subject of study in the /c^-factorization 
approach. 

Therefore, in the present paper in addition to the KMR approach we propose a some 
simplified way to evaluate the unintegrated quark densities f q ( ) within the CCFM 

dynamics. First we convolute the CCFM-evolved gluon distribution f g (x, k. T , /z 2 ) with the 
usual unregulated leading-order DGLAP splitting function P qg (z) to obtain the unintegrated 
sea quark densities. Then we add the CCFM-evolved valence quark densities which have 
been recently evaluated and applied [19] to the jet production at the LHC conditions (in the 
framework of Monte-Carlo event generator cascade [20]). Of course, in this way we only 
simulate the last gluon splitting in the full evolution cascade and do not take into account 
contribution from quarks coming from the earlier steps of the evolution. But it is not evident 
a priori, whether the last gluon splitting dominates or not. One of the goals of our study is 
to clarify this point. In order to estimate the contribution from the quarks involved in the 
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earlier steps of the evolution we use the specific properties of the KMR approach [18] which 
enables us to discriminate between the various components of the quark distributions [21, 
22]. 

We would like to point out that, in contrast with the our previous investigation [17], the 
present study is based on the off-shell matrix elements of underlying partonic subprocesses, 
where the virtualities of both incoming gluons and quark are properly taken into account. 
Numerically, we will investigate the total and differential cross sections of the inclusive and 
associated jet prompt photon photoproduction and perform a systematic comparison of our 
predictions with the available HI and ZEUS data [2-5]. Our additional goal is to study 
specific kinematical properties of the photon-jet system which are strongly related to the 
intrinsic partonic kr- 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the /c r -factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 



2 Theoretical framework 

2.1 The subprocesses under consideration 

In ep collisions at HERA prompt photons can be produced by one of three mechanisms: 
a direct production, a single resolved production and via parton-to-photon fragmentation 
processes [23] . The direct contribution to the jp — > 7 + X process is the Compton scattering 
on the quark (antiquark) 

-f(h) + q(k 2 )^-f(p 1 ) + q(p'), (1) 

where the particles four-momenta are given in parentheses. It gives the 0(al m ) order con- 
tribution to the hadronic cross section. Here a em is Sommerfeld's fine structure constant. 
The single resolved subprocesses are 

q(h)+g(k 2 ) ^^(p^ + qip'), (2) 

9(ki) + q(k 2 )^^(p 7 ) + q(p'), (3) 
q(ki) + q{k 2 ) ->-Y(pJ+g(p'). (4) 

Since the parton distributions in a photon at leading-order have a behavior proportional 
to a em ln ji 2 1 Aq CD ~ a em /a s , these subprocesses give also the 0{a 2 em ) contributions and 
therefore should be taken into account in our analysis. 

The calculation of the off-shell matrix elements (1) — (4) is a very straightforward. Here 
we would like to only mention two technical points. First, in according to the /^-factorization 
prescription [11], the summation over the incoming off-shell gluon polarizations is carried 
with Y.^t u = kyk^/k|, where k T is the gluon transverse momentum. Second, when we 
calculate the matrix element squared, the spin density matrix for all on-shell spinors is taken 
in the standard form u(p)u(p) = p + m. In the case of off-shell initial quarks the on-shell spin 
density matrix has to be replaced with a more complicated expression [24]. To evaluate it, 
we " extend" the original diagram and consider the off-shell quark line as internal line in the 
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extended diagram. The "extended" process looks like follows: the initial on-shell quark with 
four-momentum p and mass m q radiates a quantum (say, photon or gluon) and becomes an 
off-shell quark with four-momentum k. So, for the extended diagram squared we write: 



\M\ 2 ~ Sp 

where T is the rest of the original matrix element which remains unchanged. The expression 
presented between T M and 7^ now plays the role of the off-shell quark spin density matrix. 
Using the on-shell condition u(p) u(p) = p+m q and performing the Dirac algebra one obtains 
in the massless limit m q — > 0: 

\ M \ 2 ~Jj^ f " (2k 2 p - 4(p ■ k)k) %. (6) 

Now we use the Sudakov decomposition k = xp + k T and neglect the second term in the 
parentheses in (6) in the small- a; limit to arrive at 

\ M \ 2 -^r^Pr,. (7) 

(Essentially, we have neglected here the negative light-cone momentum fraction of the in- 
coming quark). The properly normalized off-shell spin density matrix is given by xp, while 
the factor 2/xk 2 has to be attributed to the quark distribution function (determining its 
leading behavior). With this normalization, we successfully recover the on-shell collinear 
limit when k is collinear with p. 



7* 



k + m q 
k 2 — ml 



u{p)u(p) 7, 



k + m q 
k 2 — m': 



(5) 



2.2 The CCFM and KMR unintegrated parton distributions 

As it was mentioned above, in the framework of /^-factorization approach one should 
consider the unintegrated gluon and quark distributions f a ( ) instead of the conven- 

tional (collinear) parton densities a(x,/z 2 ). In the KMR approximation, the unintegrated 
quark and gluon distributions are given by the expressions [18] 



fq (xy T ,^) = T q (ki,^^ 



x 



x / dz P qq (z)-q ( -, 4 ) 9 (A - z) + P q9 (z)-g ( -, k! 
J L z \z / z \z 

x 

fg(X,\i T ,fl ) =Tg()i. T ,II 



(8) 



2tt 



x / dz 



E P M {*)-Q ( 7' k T ) + PM-9 ( 7> k r ) © (A - z 



X 



X 



(9) 



where P a b(z) are the usual unregulated leading order DGLAP splitting functions, q(x,fi 2 ) 
and g(x, fj 2 ) are the conventional quark and gluon densities, T q (k^, /i 2 ) and T 9 (k^, fi 2 ) are the 
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quark and gluon Sudakov form factors, and the theta function 0(A — z) implies the angular- 
ordering constraint A = + |kr|) specifically to the last evolution step to regulate the 
soft gluon singularities [18]. 

Another the solution for the unintegrated gluon distributions have been obtained in [25] 
from the CCFM evolution equation where all input parameters have been fitted to describe 
the proton structure function F 2 (x, Q 2 ). The proposed gluon densities (namely, sets AO and 
BO) have been applied to the number of QCD processes in the framework of the Monte-Carlo 
generator cascade [20] and in our calculations [21]. 

In the present paper we will use both these distributions in our calculations. To ac- 
complish the CCFM-evolved gluon densities, one should apply the relevant unintegrated 
quark distributions. Below we will use the following way to get the f q ( ). The un- 

integrated valence quark densities f( v \ have been obtained recently [19] from the 

numerical solution of the CCFM-like equation,]. To calculate the contribution of the sea 
quarks appearing at the last step of the gluon evolution, fi 9 \x, k^, ft 2 ), we convolute the 
CCFM-evolved unintegrated gluon distribution f g ( ) with the standard leading-order 

DGLAP splitting function P qg {z): 

f^(xy T ,^ 2 ) = ^Sii J f g (x/ z ,\4^ 2 )p qg (z)dz. (io) 

x 

Note that in the region of small k 2 ^ < q 2 the scale in the strong coupling constant a s is 
kept to be fixed at go = 1 GeV. To estimate the contribution of the sea quarks coming 
from the earlier evolution steps, / w ( ), we apply the procedure based on the specific 

properties of the KMR scheme. Modifying (8) in such a way that only the first term is 
kept and the second term is omitted and keeping only the sea quark in first term of (8) we 
remove the valence and f^ 9 '(x, k^, /i 2 ) quarks from the evolution ladder. In this way only 
the fjj s \x, ky, (j, 2 ) contribution to the f q (x, ky, /z 2 ) is taken into account. 

We would like to point out that the valence quark densities from the CTEQ 6.1 set have 
been used [19] as the starting distributions to calculate the CCFM-evolved valence quark 
distributions in a proton. However, the CTEQ collaboration does not provide the quark 
and gluon distributions in a photon (which are necessary to calculate the resolved photon 
contributions), and there is no CCFM-evolved unintegrated quark densities in a photon. 
Therefore everywhere in our numerical analysis below we apply the KMR approximation 
for the unintegrated parton densities in a photon. Numerically, in (8) and (9) we have 
tested the standard GRV-94 (LO) [26] and MSTW-2008 (LO) [27] sets of collinear parton 
densities in a case of proton and the GRV-92 (LO) [28] and CJKL (LO) [29] sets in a case 
of photon. To compare the different types of evolution, we have performed the numerical 
integration of the parton densities / a (x,k^,/i 2 ) over transverse momenta k 2 ^. In Fig. 1 we 
show the obtained "effective" valence quark distributions in a proton^ as a function of x for 
different values of /i 2 , namely /i 2 = 2 GeV 2 , /i 2 = 10 GeV 2 and /i 2 = 100 GeV 2 . The solid 
lines correspond to the CCFM-evolved unintegrated (valence) w-quark and <i-quark densities. 

1 Authors are very grateful to Hannes Jung for providing us the code for the unintegrated valence quark 
distributions. 

2 The comparison of different unintegrated gluon densities to each other can be found, for example, in [14- 
16]. 
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The dashed and dash-dotted lines correspond to the relevant KMR distributions based on 
the collinear GRV-94 (LO) and MSTW-2008 (LO) sets, respectively. We have observed some 
differences in both normalization and shape between the valence quark densities calculated 
within all these approaches. Below we will study the dependence of our numerical results 
on the evolution scheme in detail. 



2.3 Cross section for the prompt photon production 

Main formulas for prompt photon photoproduction have been obtained in our previous 
paper [17]. Here we only recall some of them. Let p e and p p be the four-momenta of the 
initial electron and proton. The direct contribution (1) to the 7£> — > 7 + X process in the 
/cT-factorization approach can be written as 

a^\ 1P ^ 1 + X) =£ / ^f q (x 2 y 2T ,^)d\4^da( iq -* ig ), (11) 

q J X 2 ^VT 

where a^q — > 7g) is the hard subprocess cross section via quark or antiquark having fraction 
X2 of a initial proton longitudinal momentum, non-zero transverse momentum k 2 T (k| T = 
—k 2T 7^ 0) and azimuthal angle (j) 2 . The expression (11) can be easily rewritten in the form 

r I 1 2 2\ i -v rrn7 n 2 ^02 # 7 

xf q (x 2 ,k 2T ,n )dy 1 dE+dk 2T — — , 

where \M.{ r yq — > 7g) \ 2 is the hard matrix element squared which depends on the transverse 
momentum k| T , s = (A7 +p p ) 2 is the total energy of the subprocess under consideration, y 1 , 
Ej. and 7 are the rapidity, transverse energy and azimuthal angle of the produced photon 
in the 7£> center-of-mass frame, and a = Ej, exp y 1 / y/s. 

The formula for the resolved contribution to the prompt photon photoproduction in the 
/cT-factorization approach can be obtained by the similar way. But one should keep in mind 
that the convolution in (11) should be made also with the unintegrated parton distributions 
f2(x, k^, n 2 ) in a photon, i.e. 



x / — f b (x 2 ,kl T , /i 2 )dkl T ^p-dcr(ab -> 7c), 

J X2 ^7T 



a < b Xl (13) 



where a,b,c = q and/or g, a(ab — > 7c) is the cross section of the photon production in the 
corresponding parton-parton interaction (2) — (4). Here parton a has fraction x\ of a initial 
photon longitudinal momentum, non-zero transverse momentum kiT {k 2 T = —k\ T 7^ 0) and 
azimuthal angle <fi\. We can easily obtain the final expression from equation (13). It has the 
form 

^( 1p ^ + X) = E/^|M(«M 1c )Px 

a ,b \ > 14 

9 9 9 9 9 9 rlrh, rlrh„ rlrhH V ' 

xf2(x u k( T , n 2 )f b (x 2 , k 2 2T , /2 2 )dk 2 1T dk 2 2T dE^dy^dy 



2tt 2tt 2tt 
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where y c is the rapidity of the parton c in the 7p center-of-mass frame. It is important that 
the hard matrix elements squared \Ai(ab — > 7c) | 2 depend on the transverse momenta \t 2 T 
and k 2 T . We would like to note that if we average the expressions (12) and (14) over (f>i and 
4>2 and take the limit ~k\ T — > and — > 0, then we obtain well-known expressions for the 
prompt photon production in leading-order (LO) perturbative QCD. 

The experimental data [2-5] taken by the HI and ZEUS collaborations refer to prompt 
photon production in the ep collisions, where the electron is scattered at small angle and 
the mediating photon is almost real (Q 2 ~ 0). Therefore the 'jp cross sections (12) and (14) 
need to be weighted with the photon flux in the electron: 

da(ep -> e' + 7 + X) = J f l/e (y)da{jp ->• 7 + X)dy, (15) 

where y is a fraction of the initial electron energy taken by the photon in the laboratory 
frame, and we use the Weizacker- Williams approximation for the bremsstrahlung photon 
distribution from an electron: 

a em f l + {l- yf Q maxj _ 9 2 ( 1 1 \\ 



Here m e is the electron mass, Q min = m 2 e y 2 / (1 — y) 2 and Q max = 1 GeV 2 , which is a typical 
value for the photoproduction measurements at HERA. 

The multidimensional integration in (12), (14) and (15) has been performed by means of 
the Monte Carlo technique, using the routine VEGAS [30]. The full C++ code is available 
from the authors on request 

2.4 Fragmentation contributions and isolation 

In order to reduce the huge background from the secondary photons produced by the de- 
cays of 7T°, 1] and to mesons the isolation criterion is introduced in the experimental analyses. 
This criterion is the following. A photon is isolated if the amount of hadronic transverse 
energy E)p d , deposited inside a cone with aperture R centered around the photon direction 
in the pseudo-rapidity and azimuthal angle plane, is smaller than some value E^ ax : 

rrihad ^ nmax 

T ~ T ' (17) 

The both HI and ZEUS collaborations take R = 1, E™ ax = eEj, with e = 0.1 in the exper- 
iments [2-5]. The isolation criteria not only reduces the background but also significantly 
reduces the fragmentation components. It was shown [7, 8] that after applying the iso- 
lation cut the contribution from the fragmentation subprocesses is only about 5 or 6% of 
the total prompt photon cross section. Therefore in our further analysis we will neglect 
the small fragmentation contribution and consider only the direct and resolved production 
mechanisms. 



3 lipatov@ theory, sinp . msu.ru 
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3 Numerical results 



We now are in a position to present our numerical results. First we describe our theo- 
retical input and the kinematical conditions. After we fixed the unintegrated parton distri- 
butions in a proton and in a photon, the cross sections (12) and (14) depend on the energy 
scale /i. As it often done [7, 8] for prompt photon production, we choose the renormalization 
and factorization scales to be /i = £Ej,. In order to estimate the scale uncertainties of our 
calculations we will vary the parameter £ between 1/2 and 2 about the default value £ = 1. 
We use LO formula for the strong coupling constant a s (fi 2 ) with rij = 4 massless quark 
flavours and A QC d = 200 MeV, such that a s (Mf) = 0.1232. 

3.1 Inclusive prompt photon photoproduction 

Experimental data for the inclusive prompt photon production at HERA come from both 
the ZEUS and HI collaborations^. Two differential cross section are determined: first as a 
function of the transverse energy Ej>, and second as a function of pseudo-rapidity rp . The 
ZEUS data [2] refer to the kinematic regional defined by Ej, > 5 GeV and —0.7 < rp < 0.9 
with electron energy E e = 27.5 GeV and proton energy E p = 820 GeV. The fraction y of the 
electron energy trasferred to the photon is restricted to the range 0.2 < y < 0.9. Additionally 
the available ZEUS data for the prompt photon pseudo-rapidity distributions have been given 
also for three subdivisons of the y range, namely 0.2 < y < 0.32 (134 < W < 170 GeV), 
0.32 < y < 0.5 (170 < W < 212 GeV) and 0.5 < y < 0.9 (212 < W < 285 GeV). The more 
recent HI data [5] refer to the kinematic region defined by 5 < E^ < 10 GeV, — 1 < rp < 0.9 
and 0.2 < y < 0.7 with electron energy E e = 27.6 GeV and proton energy E p = 920 GeV. 

The transverse energy and pseudo-rapidity distributions of the inclusive prompt photon 
production for different kinematical regions are shown in Figs. 2 — 4 in comparison with 
the available HERA data [2, 5]. The solid histograms correspond to the results obtained 
using the KMR approximation for the unintegrated quark and gluon densities in a proton 
and in a photon (supplemented with the GRV-94 and GRV-92 parametrizations, respec- 
tively). The dashed and dash-dotted histograms correspond to the results obtained with 
the CCFM-evolved unintegrated quark fjj v \x, kf, /i 2 ), fjf\x, k 2 -, /j 2 ) and gluon f g ( 
distributions in a proton, as it was described in Section 2.2. These calculations are based on 
the CCFM set AO and set B0 gluon densities. The numerical predictions have been obtained 
by fixing both the factorization and normalization scales at the default value /i = Ej>. One 
can see that the HI and ZEUS data [2, 5] can be reasonably well described by using the KMR 
unintegrated parton densities. This is in a full agreement with our previous observations [17]. 
Our predictions tend even to slightly overshoot the ZEUS data at high values of y variable 
and large photon pseudo-rapidity rp (see Fig. 4). Concerning the CCFM predictions, the 
results coming from the CCFM and KMR parton densities are very similar to each other 
in the forward region, rf 1 > 0.4. However, we find the some underestimation of the HERA 
data in the rear pseudo-rapidity rp region. One of the possible reasons of such disagree- 

4 Very recently the HI collaboration has presented the data [6] which have been analysed in the kr- 
factorization approach supplemented with the KMR partons. 

5 Here and in the following all kinematic quantities are given in the laboratory frame where positive OZ 
axis direction is given by the proton beam. 
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uPDF (proton) + uPDF (photon) 


cr (HI region) [pb] 


a (ZEUS region) [pb] 


KMR (GRV-94) + KMPJ(GRV-92) 
KMR (MSTW) + KMR 7 (CJKL) 
CCFM (set AO) + KMR 7 (GRV-92) 
CCFM (set AO) + KMR 7 (CJKL) 
CCFM (set BO) + KMR 7 (GRV-92) 


45.76ltSg 
36.58!^ 
38.03l|? 5 7 
33.481^1 
33.401^1 

9Q Q7+ 4 - 12 

zy.or _ 2 .63 


47.19±tH 
37.97l|i 
41.201?$ 

36.iil?:| 

36.121?;^ 

Ol.80_ 1 92 


"reduced sea" 


7.56 


7.99 



Table 1: The total cross section of inclusive prompt photon photoproduction obtained in 
the kinematic range of the HI and ZEUS experiments. The theoretical uncertainties in the 
predictions correspond to the usual scale variations, as it was described in the text. 

ment can be connected with the contributions from the sea quarks involved in the earlier 
steps of the evolution cascade (below we will refer to these contributions as to "reduced sea" 
component). Since the "reduced sea" is not taken into account in the CCFM evolution we 
use the properties of the KMR approach to perform a rough numerical estimation of this 
contribution (see Table 1 and 2), as it was described above in Section 2.2. We found that 
the "reduced sea" component gives approximately 15% contribution to the calculated cross 
sections. However, to avoid double counting we do not sum the CCFM predictions and the 
estimated "reduced sea" contributions since part of them can be already included into the 
CCFM results (via initial parton distributions which enter to the CCFM equation). 

The total cross sections of the inclusive prompt photon production are listed in Table 1. 
To study the dependence of our results on the evolution scheme we vary the unintegrated 
parton densities both in a proton and in a photon, as it was described in Section 2.2. 
Additionally we study the effect of scale variations in the calculated cross sections. We 
found that this effect is rather large: the relative difference between results for /j, = E% 
and results for \x = Ej./2 or fi = 2Ej, is about 10%. In the kinematic region of the ZEUS 
experiment our numerical predictions obtained with the KMR parton densities are rather 
close to the ones coming from the usual (based on the collinear factorization of QCD) NLO 
calculations [7, 8]. 
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3.2 Prompt photon photoproduction in association with jet 

To calculate the semi-inclusive prompt photon production rates we apply the procedure 
which has been used previously in [17]. The produced photon is accompanied by a number 
of partons radiated in the course of the parton evolution. As it has been noted in [31], 
on the average the parton transverse momentum decreases from the hard interaction box 
towards the proton. As an approximation, we assume that the parton k' emitted in the last 
evolution step compensates the whole transverse momentum of the parton participating in 
the hard subprocess, i.e. k'y ~ — kp. All the other emitted partons are collected together 
in the proton remnant, which is assumed to carry only a negligible transverse momentum 
compared to k'y This parton gives rise to a final hadron jet with = |k'r| in addition to 
the jet produced in the hard subprocess. From these hadron jets we choose the one carrying 
the largest transverse energy, and then compute the cross section of prompt photon with an 
associated jets. 

The experimental data for this process were obtained by the HI and ZEUS collaborations. 
The HI collaboration presented the cross sections [5] measured differentially as a function of 
£y\ and the pseudo-rapidities rp and ?7 jct in the kinematic region defined by 5 < E% < 10 
GeV, E^ > 4.5 GeV, -1 < < 0.9, -1 < rj iet < 2.3 and 0.2 < y < 0.7 with electron energy 
E e = 27.6 GeV and proton energy E p = 920 GeV. The more recent ZEUS data [4] refer to 
the kinematic region defined by 5 < E% < 16 GeV, 6 < E^ < 17 GeV, -0.74 < rp < 1.1, 
— 1.6 < rf ct < 2.4 and 0.2 < y < 0.8 with the same electron and proton energies. 

The results of our calculations are shown in Figs. 5 — 8 in comparison with the HERA 
data. One can see that the situation is very similar to the inclusive production case. The 
dstributions measured by the HI collaboration are reasonably well reproduced by our calcu- 
lations supplemented with the KMR unintegrated parton densities. However, there is some 
discrepancy between the predictions and the ZEUS data. It seems that the origin of this 
disagreement is connected with the lowest bin in the Ej, distribution, where our theoretical 
results are about 2 times below the ZEUS measurements (see Fig. 5, right panel). In order to 
investigate it in more detail, we have repeted the calculations with an additional cut on the 
photon transverse energy, namely Ej. > 7 GeV (keeping the other cuts the same as before). 
Our results compared to the ZEUS data are shown in Fig. 9. We found a perfect agreement 
between the theoretical predictions (based on the KMR parton densities) and the data after 
applying this additional cut (see also [4]). Note that the KMR-based results agree with the 
HI measurements [5] in a whole Ej, range. 

Concerning the CCFM predictions, we found again that they are below the HERA data. 
In our opinion, it is connected with the missing "reduced sea" component (which gives about 
20% contribution to the total 7 + jet cross section, see Table 2). Note also that the shape 
of all predicted pseudo-rapidity r] ict distributions (based on the CCFM as well as on the 
KMR unintegrated parton densities) coincide with the ones calculated in the collinear NLO 
pQCD approximation [7, 8]. As it was pointed out [5], the shape of this distribution is not 
reproduced well by the LO pQCD calculations. This fact demonstrates that the main part 
of the collinear high-order corrections is already included at LO level in /c T -factorization 
formalism (see also [14-16] for more information). 

Now we turn to the total cross section of the prompt photon and associated jet photopro- 
duction at HERA. Results of our calculations within the framework of the fey-factorization 
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Source 


er(7 + jet) [pb] (region I) 


(7(7 + jet) [pb] (region II) 


ZEUS measurement [4] 
NLO QCD [7] 
NLO QCD [8] 

KMR (GRV-94) + KMR 7 (GRV-92) 
KMR (MSTW) + KMR 7 (CJKL) 
CCFM (set AO) + KMR 7 (GRV-92) 
CCFM (set AO) + KMR 7 (CJKL) 
CCFM (set BO) + KMR 7 (GRV-92) 

OOrlvl [SGI D\J) + rvMit, ' [vj JiyIjJ 


33.1 ±3.0 (stat.)_J;S (syst.) 
23.3±}; 9 7 
23.5_];;g 
23.1011^ 

19.28io;s9 
17.131};! 
15.29l?:g| 
15.68io!68 

i q o r +0.81 
10.00_o.82 


13.8 ± 1.2 (stat.)_J;e (syst.) 

i4.9±i;g 
i3.4i : 9 

14.881^ 

1 q+0.38 

ii.n±_:E 
io.o6±8:g 
io.26±8:fg 

y - 1U -0.40 


"reduced sea" 


4.49 


3.11 



Table 2: The total cross section of prompt photon and associated jet photoproduction ob- 
tained in the kinematic range Q 2 < 1 GeV 2 , 5 < ET T < 16 GeV, 6 < E^ < 17 GeV, 
-0.74 < rp < 1.1, -1.6 < V et < 2.4 and 0.2 < y < 0.8 (region I). An additional cut 
Ej. > 7 GeV is applied in the region II. 

approach compared to the ZEUS experimental data [4] are listed in Table 2. Similar to the 
inclusive photon production case, in these calculations we study the dependence of the pre- 
dicted cross sections on the evolution scheme and the relative effects of scale variations. The 
measured cross sections are described reasonably well using the /^-factorization approach 
and the KMR-constructed unintegrated parton densities. 

The most important variables for testing the structure of colliding proton and photon 
are the longitudinal fractional momenta of partons in these particles. In order to reconstruct 
the momentum fractions of the initial partons from measured quantities the observables x° bs 
and a;° bs are introduced in the ZEUS analysis [3, 4]: 

r o bs _ BT T e^ + E^e-^ bs E^ + E^ 

^ ~ 2yE e ' X * ~ 2E P ' [L7) 

The x° bs distribution is particularly sensitive to the photon structure function. It is known 
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that at large x° region (x° > 0.85) the cross section is dominated by the contribution 
of processes with direct initial photons, whereas at x° bs < 0.85 the resolved photon contri- 
butions dominate [4, 5]. Instead of using the x° bs and x° hs variables, the HI collaboration 
refers [5] to x^° and x^° observables given by 

lo _ EUe^+e-n _ E^+e^) 

X i ~ 2yE e ' X ? ~ 2E P " ^ 

It was argued [5] that these quantities make explicit use only of the photon energy, which is 
better measured than the jet energy. Our predictions for all these observables compared to 
the HI and ZEUS data [4, 5] are shown in Figs. 10 and 11. We conclude again that KMR 
predictions reasonable agree with the HERA data for both direct and resolved production 
mechanisms. The sizeble contribition from the "reduced sea" quarks appears only for the 
direct production and practically negligible for the resolved one. 

Further understanding of the process dynamics and in particular of the high-order correc- 
tion effects may be obtained from the transverse correlation between the produced prompt 
photon and the jet. Specially the HI and ZEUS collaborations have measured [3-5] the 
distribution on the component of the prompt photon's momentum perpendicular to the jet 
direction in the transverse plane, i.e. 

p± = \pl x p J r |/|P J ?I = #?sinA0, (19) 

where A0 is the difference in azimuth between the photon and the accompanying jet. The 
ZEUS collaboration have measured [3] also the distribution on the A<j) angle. In the collinear 
leading-order approximation, these distributions must be simply delta functions S{pj_) and 
5(0 — 7r), since the produced photon and the jet are back-to-back in the transverse plane. 
Taking into account the non-vanishing initial parton transverse momentum leads to the 
violation of this back-to-back kinematics in the /^-factorization approach. The normalised 
p± and A0 distributions compared to the HI and ZEUS data [3-5] are shown in Figs. 12 and 
13 separately for the regions x^° < 0.85 and x^° > 0.85 (in the case of ZEUS measurements 
for x° bs > 0.9 only). One can see that both the CCFM and KMR predictions are consistent 
with the data for all p± values at x^° > 0.85 (or x° hs > 0.9) and tend to underestimate the 
data in the large region at x^° < 0.85. However, this underestimation is not significant 
and therefore we can conclude that the CCFM-evolved parton densities reasonably well 
simulates the intrinsic partonic k^- The /^-factorization predictions depicted in Fig. 11 are 
very similar to the ones [7] obtained in the collinear factorization of QCD at NLO level. 
The NLO calculations performed by another group [8] give a better description of the p± 
distributions at x^° < 0.85 than the ones [7] since in this kinematical region the cross 
section is dominated by 0(a s ) corrections to the processes with resolved photons, which are 
not included in the calculations [7]. 

As a final point, we should mention that the corrections for hadronisation and multiple 
interactions have been taken into account in the NLO analysis of the available HERA data [2- 
5] performed in the framework of collinear factorization of QCD. The correction factors are 
typically 0.8 — 1.2 depending on a bin. These corrections are not taken into account in our 
consideration. 
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4 Conclusions 



In the present paper the evaluated CCFM and KMR unintegrated quark and gluon densi- 
ties have been applied to the analysis of the recent experimental data on the prompt photon 
photoproduction taken by the HI and ZEUS collaborations at HERA. Our consideration 
is based on the off-shell matrix elements of underlying partonic subprocesses (where the 
transverse momenta of both quarks and gluons are properly taken into account) and covers 
both inclusive and associated with the hadronic jet production rates. We have studied the 
dependences of our numerical results on the evolution scheme and on the standard scale 
variations. To evaluate the unintegrated quark densities within the CCFM dynamics we 
have calculated separately the contribution of valence quarks, sea quarks appearing at the 
last step of the gluon evolution and sea quarks coming from the earlier gluon splittings. In 
first time the contribution from the last gluon splitting has been calculated as a convolu- 
tion of the CCFM-evolved unintegrated gluon distribution with the standard leading-order 
DGLAP splitting function P qg (z). The contribution from the sea quarks involved into the 
earlier evolution steps has been estimated in the framework of the KMR approximation. 

We have found a reasonable agreement between our predictions and the available data. 
The contributions to the total photon cross section from the quarks emerging from the earlier 
steps of the parton evolution rather than from the last gluon splitting are estimated to be of 
15 - 20% approximately. Additionally we have studied the specific kinematical properties of 
the photon-jet system which are strongly sensitive to the transverse momentum of incoming 
partons. We have demonstrated that the /c T -factorization approach supplemented with the 
CCFM and KMR parton dynamics reasonably well simulates the intrinsic partonic k T - 

Note that in our analysis we neglect the contribution from the fragmentation processes 
and from the direct box diagram (7*7 — > 7(7). As it was claimed in [7], the direct box diagram, 
which is formally of the next-to-next-to-leading order, gives approximately 6% contribution 
to the total NLO cross section. The problem of taking into account the contribution of box 
diagrams with initial off-shell gluons in the framework of /^-factorization is still open. 
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Figure 1: The effective (k^-integrated) valence quark distributions in a proton as a function 
of x for different values of jj?. The solid lines correspond to the CCFM-evolved quark 
distributions. The dashed and dash-dotted lines correspond to the KMR predictions based 
on the collinear GRV-94 (LO) and MSTW-2008 (LO) sets, respectively. 
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Figure 2: The differential cross sections da/dE^ for the inclusive prompt photon photopro- 
duction at HERA. The solid histograms correspond to the results obtained using the KMR 
quark and gluon densities in a proton and in a photon. The dashed and dash-dotted his- 
tograms correspond to the results obtained with the CCFM-evolved quark f^ v '(x, k^, // 2 ), 
fjf\x, k|, fx 2 ) and gluon f g (x, kf,, /i 2 ) distributions in a proton. In these calculations we use 
CCFM set AO and set BO gluons, respectively. The experimental data are from HI [5] and 
ZEUS [2]. 
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Figure 3: The differential cross sections da/drf for the inclusive prompt photon photopro- 
duction at HERA. Notation of all histograms is the same as in Figure 2. The experimental 
data are from HI [5] and ZEUS [2]. 
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Figure 4: The differential cross sections da/drf for the inclusive prompt photon photopro- 
duction at HERA calculated in the different kinematical regions. Notation of all histograms 
is the same as in Figure 2. The experimental data are ZEUS [2]. 
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Figure 5: The differential cross sections da/dE^ for the prompt photon + jet production at 
HERA. Notation of all histograms is the same as in Figure 2. The experimental data are 
from HI [5] and ZEUS [4]. 
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Figure 6: The differential cross sections da/drf for the prompt photon + jet production at 
HERA. Notation of all histograms is the same as in Figure 2. The experimental data are 
from HI [5] and ZEUS [4]. 
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Figure 7: The differential cross sections da/dE? for the prompt photon + jet production 
at HERA. Notation of all histograms is the same as in Figure 2. The experimental data are 
from HI [5] and ZEUS [4]. 
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Figure 8: The differential cross sections da/d7f et for the prompt photon + jet production at 
HERA. Notation of all histograms is the same as in Figure 2. The experimental data are 
from HI [5] and ZEUS [4]. 
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Figure 9: The differential cross sections da /drf, da jdrf et and dajdE^ for the prompt photon 
+ jet production at HERA. The additional cut Ej,>7 GeV has been applied. Notation of 
all histograms is the same as in Figure 2. The experimental data are from ZEUS [4]. 
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Figure 10: The differential cross sections da/dx^° and da/dXp° for the prompt photon + jet 
production at HERA. Notation of all histograms is the same as in Figure 2. The experimental 
data are from HI [5]. 




Figure 11: The differential cross sections da/dx^ s for the prompt photon + jet production at 
HERA. Notation of all histograms is the same as in Figure 2. The additional cut E% > 7 GeV 
has been applied on the right panel. The experimental data are from ZEUS [4]. 



21 




P T (GeV) 




p T (GeV) 



Figure 12: The normalized differential cross sections 1/a da/dp± for the prompt photon 
+ jet production at HERA. Notation of all histograms is the same as in Figure 2. The 
experimental data are from HI [5]. 
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Figure 13: The normalized differential cross sections 1/a da/dp± and l/ada/dA(f) for the 
prompt photon + jet production at HERA. Notation of all histograms is the same as in 
Figure 1. The experimental data are from ZEUS [3]. 
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